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We present a Gaussian state analysis of the entanglement generation between two macroscopic 
atomic ensembles due the continuous probing of collective spin variables by optical Faraday rotation. 
The evolution of the mean values and the variances of the atomic variables is determined and 
the entanglement is characterized by the Gaussian entanglement of formation (GEoF) and the 
logarithmic negativity. The effects of induced opposite Larmor rotation of the samples and of light 
absorption and atomic decay are analyzed in detail. 



INTRODUCTION 



Macroscopic samples of atoms as a resource of entan- 
glement have attracted attention because of their robust- 
ness to single particle losses which leads to msec lifetimes 
of the entangled states and because of the effective cou- 
pling to light as needed for quantum repeaters and memo- 
ries in quantum communication networks 0, 0] ■ The the- 
oretical proposal 0,0] surprisingly showed that by merely 
probing the state of atomic samples with light from a 
classical light source, one induces an atomic dynamics 
where the quantum state evolves by state reduction to 
entangled states. The experimental implementation of 
the proposal 0] led to the first demonstration of entan- 
glement between macroscopic (~ 10 12 ) numbers of atoms 
In this work we extend the theoretical Heisenberg 
picture analysis in 0] with an analysis addressing directly 
the quantum state of the atoms and its time evolution due 
to the interaction with the continuous wave (cw) probe 
field, the back action of the measurements, obtained con- 
tinuously in time, and light absorption and atomic decay. 
We note that a quantum trajectory approach with sim- 
ulated state vector dynamics was presented in 0] to 
provide a microscopic description of the dynamics, but 
because of the dimensions of the Hilbert spaces involved, 
these simulations were restricted to a few tens of atoms. 
We retain in this work the careful attention of 0] to 
the quantum mechanical effects of the measurement in a 
treatment of macroscopic samples by a practically exact 
Gaussian Ansatz for the quantum states. This permits 
the use of the powerful formalism of correlation matrices 
for Gaussian states [flllicj. 

The successful experimental verification of entangle- 
ment of macroscopic samples utilizes Larmor rotation 
of the samples. This presents an experimental advantage 
compared to measuring on only one EPR quadrature at a 
time, but until now there has been no thorough theoreti- 
cal examination of the precise effect of these rotations on 
the entanglement generation rate. We give explicit ana- 
lytic expressions for the entanglement generation rate as 
a function of rotation frequency in the absence of light 



absorption and numerical results in the presence of light 
absorption and atomic decay. 

In section|n]we introduce the description of the atomic 
and light variables and briefly discuss the basic interac- 
tion. In section ITTT1 we introduce a Gaussian description 
of the interaction and solve a nonlinear differential equa- 
tion for the atomic variables in the absence of light ab- 
sorption. The entanglement is quantified in terms of the 
GEoF. In section IW1 we describe the effects of light ab- 
sorption and we present analytic solutions for the atomic 
variables in the case of small decoherence effects and nu- 
merical solutions for the general case. In section El we 
discuss the evolution of the mean values of the atomic 
variables during the interaction. Section IVII concludes 
the paper. 



II. SETUP AND INTERACTION 
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Figure 1: A continuous wave light beam linearly polarized 
along the x-axis is sent through two macroscopic samples of 
atoms optically pumped with collective spins in the positive 
and the negative x-direction respectively. The polarization 
rotation of the field is monitored continuously. A bias field 
along the x-direction is applied to induce Larmor rotation of 
the atomic spins in the y-z plane during the measurement. 

We consider the system studied experimentally in 0,0] 
and sketched in Fig. Qwith two macroscopic samples of 
spin 1/2 particles polarized along the positive and neg- 
ative x-axis respectively. The samples interact with an 
off resonant linearly polarized light beam giving rise to 
a Kerr-interaction between the macroscopic spin oper- 
ator, J, and the Stokes operator of light, S. We as- 
sume the two atomic samples to be prepared close to 
the maximally polarized state along x with magnitudes 
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Jxi = Jx = N a /2 and J X 2 = —J x , where N a is the com- 
mon number of atoms in each of the two samples. Here 
and in the remaining part of this paper we set fi — 1. For 
macroscopic samples of atoms the quantum mechanical 
uncertainty in J x is negligible compared to the magni- 
tude of J x which can therefore be considered a classical 
number. With large number of photons in the probing 
beam the same argument applies to the Stokes vector 
component S x . Defining a vector of observables 
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the system is described to a good approximation by oper- 
ators which obey the usual position and momentum com- 
mutation relation [x^Xj] = [pi,Pj] — 0,[xi,pj] — iSij. 
The Hamiltonian for the interaction between the light 
and either of the two samples will be given by an expres- 
sion of the form Hi — npiPL- The coupling strength k is 
proportional to the square root of the number of atoms 
and the square root of the photon number, and it will 
be related to other physical parameters in the numerical 
examples presented below. 

To model the continuous interaction between the 
atoms and the incoming cw-light field we propose along 
the lines of ^lj to split the light field into independent 
slices of duration r. The interaction between the sam- 
ples and each light segment are then treated one after 
the other. The continuous interaction and detection of 
the resulting field then corresponds to taking the r — > 
limit. 

Since both the number of photons and atoms are very 
large and the initial polarized state of the atoms and the 
light is a minimum uncertainty state, a Gaussian distri- 
bution function for the quantum variables is valid. This 
form is preserved both by the interaction and by the de- 
tection so we can use the powerful formalism of 
correlation matrices to describe the dynamical evolution 
of the system. Within the Gaussian approximation all 
information is contained in the first two moments of the 
quantum variables. We are interested in the entangle- 
ment properties of the samples which are not changed by 
local displacement operations so the second moments are 
of primary interest. These are collected in the 6x6 covari- 
ance matrix defined by j i} = 2Re {(yi - (yi)){y 3 ~ (Uj)))- 
Knowing how the covariance matrix is updated during 
the interaction and by the detection allows us to monitor 
the dynamics real-time. 



III. NO LIGHT INDUCED DECOHERENCE 

We model the evolution of the atomic system from t to 
t + t by taking sequentially into account the interaction 



of initially coherent light with each atomic sample, the 
rotation of the samples, and the homodyne detection of 
the xl quadrature of the light. The evolution of y in the 
Heisenberg picture due to the interaction between the 
light segment and samples 1(2) is given by y(t + r) = 
Si(2)y(t) with the interaction matrices: 
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Taking into account also the Larmor precession and the 
detection, the correlation matrix will evolve according to: 



7 (t + r) = M[R ■ S 2 ■ Si ■ j(t) ■ Sf ■ £T • R 1 



(2) 



where R denotes a block diagonal matrix rotating the 
atomic variables of the samples an angle ±wt and leaving 
the light variables unchanged. M[...] denotes the effect 
of the homodyne detection. Let the covariance matrix 
before the homodyne measurement be given by: 
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(3) 



where j a is a 4x4 matrix describing the atomic subsys- 
tem, and 7b is a 2x2 matrix describing the light system. 
All atom-light correlations are contained in j c . After 
the detection the atomic part of the correlation matrix is 
then given by 0, 0] : 

la -> la ~ 7c(7T7 6 7r)~7j (4) 

where it — diag(l, 0) because one quadrature of the light 
is assumed to be detected perfectly and ( • )~ denotes the 
Moore-Penrose pseudo-inverse of a matrix. 

When Eq. (J2J is evaluated for short time segments r, 
the change in 7 is quadratic in k t . k% is proportional 
to the photon number in the beam segment, i.e. propor- 
tional to t and rewriting t, the differential limit 
for the atomic correlation matrix can be formed: 

dla 

dt 

where: 



-JT = r 7a + laV T + k 2 (A- 7a S 7a T ) 



(5) 



.4 








1 































1 





i) 







1 


i 












\o 








0/ 






1 








(6) 



and 






to 








— LU 




















— LU 








UJ 






(7) 



We note that the evolution of the atomic covariance 
matrix caused by the measurements given by Eq. (@} is 
deterministic (despite the random outcomes of the detec- 
tion) and non-linear. 
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A. Ricatti equation and solution 

The nonlinear differential Eq. <j£| can be solved using 
the Ricatti method as eg. mentioned in the appendix of 
[l2T |. The generic Ricatti equation is: 



dV 
~~dt 



= C DV(i) - V(t)A - V(t)BV(t) (8) 



Using the decomposition V(t) = W(t)U _1 (i) it can 
be shown that the nonlinear differential equation can be 
replaced by the linear equation: 




D C 
B A 



{ U(t) 



(9) 



Matching our equation to the generic Ricatti equation 
and observing that 7J = "f a we obtain the linear set of 
equations 




W(i) 
U(t) 



(10) 



We see that as 9 — > we get 1 + 2/tf and 1, i.e. the 
characteristic squeezing and anti-squeezing of the mea- 
sured and their conjugate variables and no change in the 
unobserved ones. For 9 — > 00, i.e. after many rotations 
we squeeze the two quadratures symmetrically and anti- 
squeeze their conjugate variables by the factor 1 + k 2 , 
Note that the reduced squeezing by a factor of two in 
the rotated case comes from the fact that we effectively 
only spend half the time measuring on each quadrature. 
The result without rotations matches that of Q and the 
strongly rotated result agrees with calculations from 0. 



B. Gaussian Entanglement of Formation 

As an entanglement measure we choose the recently 
proposed Gaussian Entanglement of Formation(GEoF) of 
131. This measure agrees with the Von Neuman entropy 
for pure states and it can easily be calculated from the 
covariance matrix. It is given by: 

GEoF(A) = c + (A)log 2 [ c+ (A)] - c_(A)log 2 [c_(A)] 

(14) 

where c±(A) = {A- 1 / 2 ± A 1 / 2 ) 2 / '4 and A 2 = Varfci - 
X2)Var(pi +P2). Note that the small A approximation: 



where we have used that r = — r . 

Choosing the W and U matrices to start out as 4x4 
identity matrices this system of coupled linear differential 
equations can be solved. The result is fairly complicated 
but can be simplified by applying a time dependent ro- 
tation of =Fu;i/2 to sample one and two respectively. A 
further simplification can be made by noting that the 
measured quadratures are really the sum of p's and the 
difference of x'es. In the sum/difference basis: 



f;EoF(A,^lu, 2 (-) +I -L 



(15) 



shows an error of 10" 5 at A = 1/100, 0.001 at A = 1/10, 
and only 1% at A = 1/5 so it is widely applicable. From 
Eqs. and l(T3j) A 2 can be shown to be 
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where the first and the third columns are the basis vectors 
corresponding to the sum of x'es and p's and the second 
and the fourth correspond to the difference, we get the 
sum/ difference correlation matrix: 
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where a± = 1 + k 2 t ± sm(u>t). Defining the total accu- 
mulated interaction at time t t and the total 
rotated angle 9 = uit we get the main result of this sec- 
tion: 



a± = 1 + k 2 ± -± sin(6>) 
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which tells us that the entanglement will scale as 
log 2 (K 2 ) = 21og 2 (/i t ) with rotation and with log 2 (/Ct) 
without rotations. The factor of two is exactly what one 
should expect since rotations enable us to squeeze both 
quadratures compared to squeezing only one of them. 
This effect was also observed in p|. 

In Fig. [21 we show the entanglement plotted as a func- 
tion of Kt and 9. As can be seen, the transition from the 
static to the rotated regime occurs before one full revolu- 
tion is reached. That is, given a total interaction time, t, 
there is no real gain in choosing the frequency larger than 
^crit = 1/t. This arises from the fact that we measure 
x- cos(9) + p + sin(0) and this operator commutes with 
the operators measured at all previous times (other val- 
ues of 9). It is therefore the accumulated measurement 
on each quadrature that counts and this will not benefit 
from many rotations compared to a single rotation. In 
the experiments of 0, 0] vl ~ 320kHz and T pro b c ~ 1ms 
so the spins will rotate several hundreds of times in one 
pulse and the results are thus firmly obtained in the lu 
independent regime. Note that both axes in Fig. rep- 
resent time dependent quantities. 
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Figure 2: 3D plot of GEoF as a function of the accumulated 
interaction strength Kt and the rotated angle 0. The plot 
clearly shows that the transition from the static to the rotated 
regime has occurred already before one full revolution of the 
atomic spins has taken place. In both the static and the 
rotated regions we clearly see the logarithmic behavior of the 
GEoF as a function of Kt. 



IV. LIGHT ABSORPTION AND ATOMIC 
DECAY 

Due to the interaction with the light field every atom 
has a rate by which it is excited and decays by sponta- 
neous emission to any of the atomic ground states, and 
every photon is absorbed with a given probability. This 
atomic depumping parameter and the photon absorp- 
tion probability are given by r] T = t\t — 3>t-§ (■£) , and 

e = N a ^ (^) 2 respectively. $ is the photonic flux, A is 
the cross section of the atomic sample illuminated by the 
light, A is the detuning from resonance, a is the cross sec- 
tion on resonance for the probed transition, and T is the 
corresponding spontaneous decay rate. The optical den- 
sity on resonance is ao = N a ^. This gives the relation for 
the coupling parameter introduced above, k^(0) = r] T ao. 
The decay of the mean spin makes the coupling constant 
time dependent, (0)e vt . We will in our numer- 

ical simulations use the experimentally motivated values 
of 5 MHz for the decay rate T and 1000 MHz for the de- 
tuning A. Many results will be presented as a function 
of ao so a brief discussion of the experimentally realiz- 
able optical densities is in order. In the experiments of 
ao ~ 5 whereas optical densities of 100 or more can 
be routinely achieved in MOTs (see e.g. 0]). Using a 
Bose-Einstein condensate ao ~ 1000 can be achieved. 



A. Covariance matrix update and Ricatti equation 

We now derive an expression for the time evolution 
of the covariance matrix due to photon loss and atomic 
decay. Due to atomic decay, during a time interval r, a 



fraction r\ T of the atoms decays into a random mixture 
of the ground states, giving rise to the new value of the 
variance of one of the atomic collective spin component: 



(1 
(1 



Vt) 2 < Jl > 



-(1 - (1 - r h ) 2 )N/A 
-r) T N/4+r) r N/4 (17) 



in the i] T <C 1 limit. The atomic decay leads to a corre- 
sponding reduction of the mean spin J' x = (1 — r] T ) J x . 

At this stage the gas contains two components: the 
atoms which have not decayed, described by the first two 
terms in the latter expression and the ones which have 
decayed, described by the last term. If nothing else hap- 
pens to the atoms, subsequent interaction with the light 
can have no further effect on the random component, but 
a new fraction of atoms will be randomized and we obtain 
the iterated expression for the spin variance 



2 ,2 2 N N 
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In the limit n — ► 00 only the last term, representing the 
fully random component of the gas, will contribute. The 
geometrical series gives l/rj T , so we see that the initial 
squeezing will decay exponentially as expected and we 
will end up with N unpolarized atoms each contributing 
1 1 A to the variance as expected. Note that in the spin 1 jl 
case this coincides with the noise of atoms in the coherent 
spin state in which the atomic samples are initialized in 
this paper. This does not, however, hold for higher spin. 

Taking into account that the covariance matrix deals 
with the transverse spin components scaled by the macro- 
scopic longitudinal mean spin, we obtain for the corre- 
sponding diagonal covariance matrix element: 



la.ii = (1 - Vt) 



1 a . j 



Vt 
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(19) 

This analysis treats the atoms that have decayed and 
the ones that have not decayed on unequal footing, and 
it hence breaks with the Gaussian state Ansatz, which 
assumes that all information is in the collective variance 
and mean values for the entire atomic ensemble. The 
analysis does not make it easy to treat the coherent part 
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of the interactions and the measurement back action, 
that we are interested in, and we hence wish to inves- 
tigate, whether restoration to the Gaussian state Ansatz 
will yield a large discrepancy with the exact results. To 
this end, we go back to the update formula, Eq. Ijl7ll . and 
insist that the variance obtained here should be treated 
as the variance describing a Gaussian state ensemble, i.e., 
we do not discriminate between the two kinds of atoms. 
In subsequent time steps, we thus simply iterate the same 
expression, as if all atoms contribute evenly to the joint 
variance. The result of this iteration is readily obtained: 



la{t + T) = (1 - 7fr-)7„(t) + 



N a 



AJ x {t 



(20) 



Note that since J x (t) = (N a /2)e-^ the last term will 
diverge exponentially in time. In order to compare with 
Eq. I|19J1 we iterate the diagonal elements of Eq. 1201) n 
times: 

n 

llu = (1 - VrTllu + VrJ^i 1 ~ VtT-V. (21) 

In the continuous limit t = nr with rjr — > and n — > oo 
both Eq. iffi))) and Eq. {2H yield the same behavior: 

Ja,ii (t) = e-^a (t = 0) + sinh(77<) (22) 



This supports the use of the update formula, Eq. J2£ 
with the underlying assumption of a Gaussian state, to- 
gether with the evolution of 7 due to interaction and 
measurements. This approach was also utilized in [TlLIT^ | 
albeit not with the careful justification presented above. 

When light absorption is included the interaction of a 
light segment with one sample i is described by fl5T |: 

7(t + r) = D(r ]T ,e)S l (K T ) 1 (t)S t (K T ) T D(r lT ,e) 

+ A(f?r,e)7noiso,i (23) 

where, following the above argument, for the first sam- 
ple we have, £>i( r?i, ex) = d iag(r?x, rg, 0, 0, ex, ei), 

Di(r)i, e i) = \A ~ D i(Vi^i) and 7 no i S o,i = 
diag(£, £, 0, 0, 1, 1) (and similarly D 2 (r]2,e2) = 
diag(0,0, r?2,??2,e2,e2), 7noiso,2 = diag(0, 0, £, £, 1, 1)). 
The factor £ = N at / (J x (t)) starts out as 2 and increases 
exponentially because of the decay of the mean spin due 
to excitation and subsequent decay of atoms. 

As in section ITTTI the differential equation for the cor- 
relation matrix can be found: 



dt 



?7 + 7? T + A — (1 — e)k" i "fB"f 



,-.2. 



(24) 



where 



A = 



(k 2 +^{t)ri k 2 y/T^l \ 

£(i)77 

k 2 y/T^i K 2 +^(t)?7 

V t(t)vj 



(25) 



B = 



/ ^ 

1-e 



\o VT~^ 1 J 



(26) 



and r = r — (ry/2)I 4 , where r is defined in Eq. (0. Note 
that the coefficients are now time dependent which com- 
plicates matters slightly. 

This Ricatti equation can of course easily be solved 
numerically but is in its most general form too compli- 
cated to admit analytical solution. If noise terms arising 
from the absorption of light, i.e. all terms involving e 
are neglected, the Ricatti equation can be solved with- 
out rotations and in the strongly rotated regime. These 
solutions will be derived first and then the general results 
will be discussed. 



B. Analytical results 

1. Without rotations 

Without rotations j a becomes diagonal in the 
sum/ difference basis if we neglect the photon absorption 
(e = 0). The two components involved in the measure- 
ment give: 



T11 = Var(a; A i 
lit'- 



XA2) 

Var(pAi +PA2) 

e vt (8 cosh(S rjt) + smh(6rit)) 
Scosh(Srit) + (1 + 2a Q ) smh(5r]t) 



2a oV t} (27) 



(28) 



where S = \/l + 4«o- The two remaining components 
increase exponentially: 



722 = Var(xAi - x A 2) = lit = Var(pAi - PA2) = e vt 

(29) 

The validity of this solution can be tested by comparing 
with full numerical solutions. In Fig. 0we plot the GEoF 
as a function of k 2 (proportional with the total number 
of transmitted photons) and we see that in the case of 
dissipation, the entanglement reaches a maximum, and 
hereafter it decays and vanishes at a point when - because 
of the decay of the macroscopic spin - the atomic samples 
are in a mixed state with too few correlations to display 
actual entanglement. 

As can be seen in Fig. 03 the analytic expression fits 
very well with the nontrivial result of the numerical cal- 
culation. The comparison is made for an optical depth, 
«o = 100, corresponding to e = ■ ao — 0.0025 <^ 1. 
For larger optical depths, the quality of the analytical ex- 
pression deteriorates whereas for lower optical densities 
the error of the analytical expression is negligible. This 
thus establishes the regime of validity for the analytic 
expression. 

It is interesting to investigate when the maximum en- 
tanglement is reached and at what level as a function of 
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Figure 3: The generated entanglement vs. time of numerical 
simulation (full) and analytical solution(dashed). The optical 
density is 100 giving rise to only a slight difference between 
the two curves. The numerical solution is also represented in 
terms of the logarithmic negativity (dash-dot) which will be 
discussed in a later section. For further comparison a numer- 
ical simulation of GEoF with rotations (0totai = 100 • 2-k) is 
included (dotted). 



o;o and rj. By differentiating A 
to time we get: 



V722733 with respect 



To illustrate the use of this new description first ne- 
glect all decoherence. With the usual update equation 
for the system we can form the differential equation for 
the atomic correlation matrix: 



dt 



la 



(33) 



When this is inserted into the Ricatti equation we easily 
derive the solution: 



7«(*) = 



1 







(34) 



which coincides with the rotated solution found in (13). 

In order to treat decoherence analytically an approxi- 
mation similar to that of section lTV B 11 i.e. ignoring all e 
dependent noise contributions, has to be made. The loss 
of light in the first sample, represented by e will create 
asymmetries between the two samples which would cause 
the xa, Pa formalism of Eqs. 131132)1 to break down. The 
regime of validity is thus expected to be similar to that 
of section HVBTl i.e. a Q < 100. 

We will consider the case of small values of r)t, i.e. the 
atomic samples retain a constant polarization along the 
x-axis and the differential equation has constant coeffi- 
cients. 



arccosh( | 



-2a 2 -4a 3 +(l+5a+4 a 2 )\ 
^ 



(30) 



Note the simple inverse scaling with rj. Since rj and t 
only appear as a product in Eqs. 1(27)1 and ll29)) the inverse 
scaling of i C rit with rj means that the maximum achievable 
entanglement will be independent of rj. This behavior 
could also be predicted from the fact that as seen in Eqs. 
il28)l and l)29)) A is a function of the product r\t. 



2. Many rotations 

To solve Eq. l)24T) in the case of many rotations we 
define a new set of canonical operators 0|: 



xa 

PA 



xl = \j Jq S y (t) cos(uj L t)dt 
Jo S z (i) cos(w L t)dt 



PL 



(31) 



where J' k refers to rotating frame coordinates, i.e. coor- 
dinates rotated an angle compared to the usual lab 
frame coordinates. With these collective operators we 
regain the same formal description of the interaction as 
in the non-rotating case: 



"A 

~out —in 

l L — X L 



Pa — Pa 
Pl — Pl 



(32) 



Note that in this formalism the atomic correlation matrix 
is reduced to a 2x2 matrix whereas the correlation matrix 
of light remains a 2x2 matrix. 



la 



e-^ [(a + 4)(e" t -l) + l] 




where j3 = ^/T+TGoq. 

An expansion in small rjt yields: 



7n 



722 



(1 



-r,t){l + 
2 + rjt(7 



7(e T ' t ^-l) + (e'"' 3 + l),3 

(35) 

(36) 
(37) 



(a + i)rjt) 
+- VI + 16a ) 



2 + Tft(l + 2a + VI + 16a ) 

Eq. 13 7|l shows that a value of 1 + 2ao > 7, is required 
in order to get squeezing in pa- Again we see that the 
variances are only functions of the two variables ao and 
r]t, and in the limit of vanishing noise we obtain 711 — > 
1 + and 722 — > 1/(1 + Kf)i m accordance with our 
earlier results. 



C. General results 

We now turn to general numerical results beyond the 
regime of validity of the analytical solutions. First how- 
ever we need to note that the GEoF is not applicable in 
the presence of the large asymmetries between the two 
cells that occurs at high a n We therefore introduce the 
logarithmic negativity of [lfl |. 



1. Logarithmic negativity 

In this measure the symplectic spectrum of the partial 
transpose of the covariance matrix r ) TA is used to cal- 
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With rotations 



no rotations 
with rotations 




Figure 4: The maximum generated degree of entanglement as 
a function of ceo with and without rotations. 

culate the degree of entanglement. In the partial trans- 
position all covariance matrix elements involving pa and 
other observables are multiplied by (-1). The symplectic 
eigenvalues are most conveniently calculated by comput- 
ing the eigenvalues of the matrix ct _1 7, where a is the 
matrix specifying the commutators: o a p — — i[y a ,y/3}. 
For m systems this gives 2m complex eigenvalues A^. 
The logarithmic negativity entanglement measure is then 
given by: 

2 m 

log.ncg. = -^log 2 [min(l,2|A fe |)] (38) 
fc=i 

This is a valid entanglement monotone even for asym- 
metric samples but it does not coincide with the Von 
Neuman entropy for pure states as the GEoF does. If 
7 is diagonal in the sum/ difference basis the eigenvalues 
can be expressed in terms of the diagonal entries: 

A = ±i] jj§M + ^7 3 1 ± iiltiit - lihil) (39) 

Without losses the covariances are specified by Eq. 
(1 1 -3 1 giving a logarithmic negativity of log 2 (/-Ct ) + 1 and 
21og 2 (Kj) in the non-rotated and the rotated regimes re- 
spectively at /tf ^> 1, This is about a factor of two larger 
than the GEoF as specified in Eq. 11 "ill . The difference 
between GEoF and the logarithmic negativity in the pres- 
ence of decoherence is illustrated in Fig. where log. 
neg. (dash-dotted) is seen to be approximately a factor of 
three larger than the GEoF (full drawn). 

2. Numerical results 

In Fig. Q] we see the degree of entanglement generated 
with and without rotations. For low optical densities the 
improvement in the degree of entanglement due to rota- 
tions of the samples is close to the factor of two, which we 
got in the absence of decoherence. This is also illustrated 
in Fig. 0for the case of a = 100. For higher optical den- 
sities the benefits of rotations seem to decrease. This is 



Figure 5: The r\t at which the entanglement has decayed to 
zero with and without rotations. The solid curve is calculated 
from Eqs. Q27I291 assuming no rotations and no asymmetry 
between the two samples. 

probably the result of the complicated interplay between 
the beneficial rotations and the increasing asymmetry be- 
tween the two samples at high optical densities. It would 
certainly be a of great interest to study this further and 
to consider protocols where the samples are illuminated 
with lasers from both sides to restore their symmetry. 

In Fig. 0we show at which accumulated r\t the entan- 
glement has decreased to zero. The theoretical curve is 
calculated from Eqs. <GZJl to pflf . This solution fits the 
numerical results perfectly within the regime of validity, 
that is up to a = 200. After this point entanglement de- 
cays rapidly, indicating that light losses which introduce 
asymmetry between the two samples are detrimental to 
the entanglement. As can be seen, the entanglement is 
not only increased by rotations. The rotated samples also 
seem to be a lot less vulnerable to light losses. Note that 
the decay in all cases happens in the vicinity of rjt = 1 
which could be expected from the exponential character 
of the decay of the mean spin and the resulting exponen- 
tial growth Of 7noisc- 

V. MEAN VALUES 

The generated entanglement is of no practical use un- 
less the mean values of the atomic parameters are known. 
These will be affected by the interaction according to: 

(y(t + T))=D 2 S 2 D 1 S 1 (y(t)) (40) 

The mean values of the atomic variables after the inter- 
action (y a (t + t)) will subsequently transform according 
to the result of the homodyne measurement as: 

(y«(* + r)) -> (y„(t)) + 7c (7T7^r 1 (x, 0) T (41) 

where \ ls the difference between the measurement re- 
sult and the mean value of the detected field component. 
That is, the evolution of the mean values is determined 
by the value of the correlation matrix at a given time. 
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Note also that it is the deviation from the mean value 
of the light that determines the evolution of the atomic 
mean values and not the measurement result itself. This 
will be discussed further below. The correlation matrix 
can be made diagonal in the sum/difference basis. We 
can therefore without loss of generality assume J a {t) to 
be diagonal. 

Having established the relatively large domain of va- 
lidity of the analytical solution of the lossy system with- 
out rotations and negligible absorption losses in section 
IIV B II we now apply the results of that section to Eq. 
lUDl . Given the fact that all mean values start out as zero 
we obtain the transformation of the atomic variables: 



(y Q (* + T)) = (l-ry T /2) (y (i)) + 



X (42) 



where 733 has the value specified in Eq. II28H . We see that 
the two x-variables will have zero mean at all times and 
the two p-variables will experience the same evolution. 
After the interaction we get in the limit of vanishing e: 



755 = 2Var(ir L ) = 1 + 2 7: 



<d^2 
33 k t 



(43) 



Since we always work in the regime kJ«1 and 7I3 « 1 
in the relevant regime we can safely assume 755 = 1. 
This means that the difference between the result of the 
homodyne measurement of xl and its expectation value, 
will be a Gaussian random variable, x with zero mean 
and variance 1/2. 

Eq. 1421 can be transformed into a stochastic differen- 
tial equation with the solution: 



-y(T-t)/2 K sd 



V2 



7 3 W 



(44) 



Since (p!(T)) = (p 2 (T)), Var«pi(T)) + (p 2 (T)» will be 
four times the result of Eq. lUfill . It follows that at any 
time 



1 



73 1 + 4Var((p)) = TT ^ + 



2k? 



1 



2k? 



= 1 



(47) 



Initially, the expectation value is well determined (0) 
whereas the quantum deviation from this value is given 
by the initial Gaussian distribution. After a significant 
interaction time the quantum mechanical uncertainty will 
be reduced but the value within the initial distribution 
at which the expectation value settles is uncertain. 

Eqs. H42I44JI express the conditional mean value of the 
atomic variable in terms of the difference between the 
optical read-out and its expectation value. It is interest- 
ing to obtain similar expressions in terms of the actual 
read-out. Since the coherent light initially has zero mean 
Eq. 14011 shows us that the measured light component 
will have the mean value: 



( XL (t)) = Krdpxit)) + (p 2 (t))) 



(48) 



in the absence of e decay. We thus start by writing x m 
Eq. 1(4*2)1 as x ~ 2k t (p(t)). The atomic p variable thus 
changes as 

(p(t + T)) = (1 - V r/2) (p(t)) + K Tl S 3 Km - ^ T (p(t))) 

(49) 

where x ls the random detector output. Taking the limit 
of inifitesimal r = dt a nd defining the measured Wiener 
increment dW = \ // 2dtx, we can integrate Eq. H49|) : 

(p(T)) = [ T e -V(T-t)/2-f t T 2~ K * 7 lt(t')dt> Z**(t)dW 

Jo v2 

(50) 

In the limit of 77 = 0, 7ff(i) = 1/(1 + 2kH), and we 
can explicitly integrate the argument of the exponential 
function in Eq. lloOt . and we obtain 



where dW — \/2dtx ls a stochastic Wiener increment 
with zero mean and variance dt. We recognize an expo- 
nential memory decay for early detection events, with 
tmem ~ 2/77, and a covariance matrix weight factor, 
733 (t), based on the state of the atoms at the particu- 
lar time. From Eq 114411 it follows that the conditional 
mean value of p will be a stochastic variable with zero 
mean and variance: 



Var«p(T))) = — 



e-vir-t)^^ (45) 



To illustrate the physical interpretation of Eq. Illoll we 
neglect decoherence and rotations. In this simple case we 
can perform the integration in Eq. ll4*5|l analytically: 



Var«p(T))) 



-dt 



1 k T 



2 1 + 2k 2 T 



(46) 



(p( T )) = / 1 + 2k J JL^gdW 

i-T 

dW (51) 



H 



y/2(l + 2k 2 T) Jo 

This remarkable result shows that all measurements 
should be weighted equally in the absence of decoher- 
ence and rotations. When transformed back into regu- 
lar angular momenta the common weight factor has a 
clear interpretation as the ratio between the shot to shot 
atomic noise contribution, the so-called projection noise, 
and the total noise in complete accordance with the result 
obtained in 0. 

We stress that the result "5T|l was obtained for the non- 
rotated and non-decaying atomic systems. Decay can be 
included easily according to Eq. Q50J1 . In addition to the 
exponential damping term, this will involve a more com- 
plicated expression for 733(4), and the conditioned mean 
value of p(T) will no longer be given by the integrated 
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measurement outcome. We can understand that atomic 
mean values attained early during the measurements de- 
cay and hence they contribute less to the final (p(T)) 
than the most recent contributions to the optical detec- 
tion. Since decay is inevitable, our analysis suggests that 
experiments must be carried out so that the optical signal 
is recorded in time bins which are much shorter than the 
atomic decoherence time rf~ 1 . The results will then be 
in good agreement with our continuous update theory. 



VI. CONCLUSION 

In conclusion, we have presented a theory for the 
preparation of entangled atomic ensembles by detection 
of the Faraday polarization rotation of a continuous op- 
tical field passing though both ensembles. Our Gaussian 
Ansatz is very well justified for large atomic ensembles 
and for free space atom-light interaction, where only the 
interaction with many photons appreciably modifies the 
atomic state. The theory incorporates the interaction be- 
tween the atoms and the optical field, atomic decay, and 
the measurement induced transformation of the atomic 
state. The reduction of the full quantum state descrip- 
tion to a simple Gaussian state fully represented by a set 
of mean values and a covariance matrix makes the system 



straight forward to deal with numerically, and analytical 
results can be obtained in several important cases. 

The entanglement between the atomic ensembles is 
quantified by the Gaussian Entanglement of Formation 
and the Logarithmic Negativity, and we identify the opti- 
mal performance of the entanglement scheme in the pres- 
ence of atomic decay. Our analysis confirms a number 
of results, derived in less complete or purely numerical 
studies, and it presents an intuitive physical picture of 
the continuous transformation of the atomic quantum 
state from an initial state with no atomic correlations 
into a state with stronger correlations of the quantum 
observables, around mean values with a broader random 
distribution - but known by the experimentalist in every 
implementation of the experiment. 

The results of our analysis are relevant for current ex- 
perimental efforts to exploit entangled atomic ensembles 
for quantum purposes, but we also wish to emphasize 
the strengths of the theoretical method, which make it 
readily adapt to a wide variety of experiments. 
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